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Abstract 

We present the method of group foliation for constructing non-invariant 
solutions of partial differential equations on an important example of the 
Boyer-Finley equation from the theory of gravitational instantons. We show 
that the commutativity constraint for a pair of invariant differential opera- 
tors leads to a set of non-invariant solutions of this equation. In the second 
part of the paper we demonstrate how the hodograph transformation of the 
ultra-hyperbolic version of Boyer-Finley equation in an obvious way leads to 
its non-invariant solution obtained recently by Manas and Alonso. Due to ex- 
tra symmetries, this solution is conditionally invariant, unlike non-invariant 
solutions obtained previously. We make the hodograph transformation of 
the group foliation structure and derive three invariant relations valid for 
the hodograph solution, additional to resolving equations, in an attempt to 
obtain the orbit of this solution. 
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1 Introduction 

The standard method for obtaining exact solutions of partial differential 
equations (PDEs) by the symmetry group analysis is the symmetry reduction 
which gives only invariant solutions, i.e. the solutions which are invariant 
with respect to some subgroup of symmetry group of the PDEs PQI21E] • The 
drawback of the standard symmetry analysis is that with this method one 
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loses all non-invariant and many invariant solutions, the latter being non- 
invariant solutions of the equations obtained by symmetry reduction. 

Meanwhile, non-invariant solutions of PDEs might be very important for 
applications in physics. The classical example is the so-called K3 instanton in 
the Einstein's theory of gravitation which corresponds to the Kummer surface 
in differential geometry. The corresponding Kahler metric has no Killing 
vectors and it could be constructed once one finds non-invariant solutions 
of the complex Monge-Ampere equation, the problem which has not been 
solved for over 150 years. 

Recently we discovered that the method of group foliation can be made 
an appropriate tool for obtaining non-invariant solutions of non-linear PDEs 
that admit an infinite-dimensional symmetry group. The method involves 
the study of compatibility of given equations with a differential constraint, 
which is automorphic under an infinite dimensional symmetry subgroup, the 
latter acting transitively on the submanifold of the common solutions. By 
studying compatibility conditions of this automorphic system, i.e. the re- 
solving equations, one can provide an explicit foliation of the entire solution 
manifold of the considered equations into separate orbits. Unlike the stan- 
dard symmetry method we do not loose any solutions in the process and we 
can select orbits of non-invariant solutions by obtaining appropriate solutions 
of the resolving system. 

The idea of the method, belonging to Lie and Vessiot |%H5] , is more than a 
hundred years old being resurrected in a modern form by Ovsiannikov more 
than 30 years ago (see P and references therein). We have modified the 
method by introducing three important new ideas PITIISII^ITU] : 

1. The use of invariant cross-differentiation, involving operators of invari- 
ant differentiation and their commutator algebra for derivation of resolving 
equations and for obtaining their particular solutions. 

2. Commutator representation of resolving system in terms of the commuta- 
tors of operators of invariant differentiation. 

3. Method of invariant integration for solving automorphic system. 

We illustrate all this on a physically important example of the Boyer- 
Finley equation ^Tj which is known also as the heavenly equation 

u Z z = K(e u ) u u xx + u yy = K(e u ) u , k = ±1 (1.1) 

where u = u(t, z, z). This equation is a continuous version of the Toda lattice 
or SU (oo) Toda field [12J. It appears in the theory of gravitational instan- 
tons [T31ll4t lTo] where it describes self-dual Einstein spaces with Euclidean 
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signature having one rotational Killing vector. It also appears in the general 
theory of relativity ^H] and other physical theories. 

In the first part of this paper, on the example of the Boyer-Finley equa- 
tion, we clarify main concepts of the method including these three ideas and 
consider in detail the main steps which should be performed for obtaining 
non-invariant solutions. This approach enables us to present explicitly a fam- 
ily of non-invariant solutions of the Boyer-Finley equation depending on two 
functional parameters. 

In the second part of the paper we discuss a new non-invariant solution 
of the ultra-hyperbolic version of the Boyer-Finley equation with k = 1 
obtained recently by Manas and Alonso by the hodograph transformation of 
associated systems of the hydrodynamic type [T7] (see also [TH] for a gener- 
alization of the approach of JZI)- We point out that if one makes the hodo- 
graph transformation of the Boyer-Finley equation itself, then it becomes 
obvious how to impose a differential constraint immediately leading to the 
solution of Manas and Alonso without any use of hydrodynamical systems or 
the general theorem presented in [T7]. We show that this solution possesses 
extra symmetries depending on eight arbitrary functions of u which are not 
the symmetries of the Boyer-Finley equation. Hence, this solution is condi- 
tionally invariant ^H] though not an invariant solution. A similar argument 
shows that symmetries of the solutions obtained by Calderbank and Tod 14 
and by Martina, Sheftel and Winternitz jS] coincide with the symmetries of 
the Boyer-Finley equation itself and hence these solutions, having no extra 
symmetries, are not conditionally invariant. Finally, we make the hodograph 
transformation of the whole group foliation structure and obtain the relations 
between differential invariants valid for the solution of [17 j which together 
with the resolving equations partially describe the orbit of this solution. We 
have not yet succeeded to obtain a solution of the resolving system which 
corresponds to the orbit of the hodograph solution of [T7]. Since the addi- 
tional relations which we have found do not form a complete set of invariant 
differential constraints which uniquely select the hodograph solution, the cor- 
responding solution of the resolving system would give us the orbits of more 
general solutions of with k = 1 and also solutions with k = — 1. This 
work is still in progress. 
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2 Symmetry algebra, differential invariants 
and automorphic system 

We start with the symmetry algebra of the generators of point transformations 
for the Boyer-Finley equation [8 

T = d t , G = td t + 2d u 
X a = a(z)d z + a(z)d- z - (a'(z) + a\z))d u (2.1) 

where a(z) is an arbitrary holomorphic functions of z and prime denotes 
derivative with respect to argument. 

The Lie algebra of symmetry generators (j2.1j) is determined by the com- 
mutation relations 

[T, G] = T, [T, X a ] = 0, [G, X a ] = 0, [X a , X b ] = X aV _ ba , (2.2) 

They show that the generators X a of conformal transformations form a sub- 
algebra of Lie algebra (|2.2j) . This subalgebra is infinite dimensional since 
the generators X a depend on a(z). The corresponding finite transformations 
form an infinite dimensional symmetry subgroup of the equation (jl.lj) since 
instead of a group parameter they also involve an arbitrary holomorphic func- 
tion of z. We choose this infinite dimensional conformal group for the group 
foliation. 

Next we find differential invariants of the symmetry subgroup of con- 
formal transformations. They are invariants of all the generators X a in the 
prolongation spaces, so that they can depend on independent variables, un- 
knowns and also on partial derivatives of unknowns allowed by the prolon- 
gation order. The order of the differential invariant is defined as the order 
of the highest derivative which this invariant depends on. The highest or- 
der N of differential invariants should be larger or equal to the order of the 
equation (N > 2) and there should be n functionally independent invariants 
with n > p + q where p and q are the number of independent and dependent 
variables, respectively. Thus, we should have p = 3, q — 1, n > 4, N > 2. 

The routine calculation for N = 2 gives n = 5 functionally independent 
differential invariants up to the second order inclusively 

t, u t , u tt , p = e~ u u z2 , r, = e~ u u zt u- zt (2.3) 

and all of them are real. This allows us to express the Boyer-Finley equation 
(jl.lj) solely in terms of the differential invariants u u = up — u\. 
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The next step is choose the general form of the automorphic system. 
We choose the invariants t, u t , p as new invariant independent variables, the 
same number p = 3 as in the original equation and require that the 

remaining invariants should be functions of the chosen ones. This provides 
us with the general form of the automorphic system that also contains the 
studied equation expressed in terms of invariants ()2.3|) 

u a = Kp-v>t, r] = F(t, Ut, p) (2.4) 

The real function F in the right-hand side should be determined from the 
resolving equations which are compatibility conditions of the system ([2.4)1 . 
Then the system (12.4)) will have automorphic property, i. e. any of its solutions 
can be obtained from any other solution by an appropriate transformation 
of the conformal group. 

3 Operators of invariant differentiation and 
the basis of differential invariants 

Our next task is to find operators of invariant differentiation. They are linear 
combinations of total derivatives D t , D z , D z with respect to independent 
variables t, z, z with the coefficients depending on local coordinates of the 
prolongation space which are defined by the special property that, acting on 
any (differential) invariant, they map it again into a differential invariant. 
As a consequence, these differential operators commute with any infinitely 

oo 

prolonged generator X a of the conformal symmetry group. Being first order 
differential operators, they raise the order of a differential invariant by one. 
The total number of independent operators of invariant differentiation is 
obviously equal to 3, the number of independent variables. 

We obtain the basis for the operators of invariant differentiation PlI^ITU] 

S = D t , A = e- u u- zt D z , A = e- u u zt D- z (3.1) 

The next step is to find the basis of differential invariants which is defined 
as a minimal finite set of (differential) invariants of a symmetry group from 
which any other differential invariant of this group can be obtained by a 
finite number of invariant differentiations and operations of taking composite 
functions. The proof of existence and finiteness of the basis was given by 
Tresse j20] and in a modern form by Ovsiannikov PQ. 
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In our example the basis of differential invariants is formed by the set of 
three invariants t,u t ,p, while two other invariants u tt and rj defined in (|2.3|) 
are given by the relations 

u u = 8(u t ), 7] = e~ u u zt u 2t = A(u t ) = A(u t ) (3.2) 

All other functionally independent invariants of higher order can be obtained 
by acting with operators of invariant differentiation on the basis {t,Ut,p}. 
In particular, the following third order invariants generated from the second 
order invariant p by invariant differentiations will be used in our construction 

a = A(p), a = A(p), r = 8{p) = D t (p) (3.3) 



4 Commutator algebra of invariant 
differential operators 

The operators 8, A and A defined by (J3.1j) form the commutator algebra 
which is a Lie algebra over the field of invariants of the conformal group P . 

This algebra is simplified by introducing two new operators of invariant 
differentiation Y and Y instead of A and A and two new variables A and A 
instead of a and a, defined by 

A = rjY, A = t]Y, a = r]X, a = r]X (4.1) 

The resulting algebra has the form 

[5, Y] = (kX - 3u t - ^\ y, [5, Y] = (kX - 3u t - Y 

{Y,Y]= {T + Utp) (Y-Y) (4.2) 

With the use of operators 5, Y and Y the general form ()2.4j) of the 
automorphic system becomes 

5{u t ) =np- u 2 t , Y(u t ) = 1 (Y(u t ) = 1) (4.3) 

where the first equation is the Boyer-Finley equation and the second equation 
follows from the second relation (|3.2jl . Here we put rj = F in the equations 
(j4.1j) and in the commutation relations (|4.2j) according to the second equation 
in (|2.4j) . Then we obtain 

y = (l/F)A, Y=(1/F)A, Y{p)=X, Y(p) = X (4.4) 
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5 Resolving equations 



The system of resolving equations is a set of compatibility conditions between 
the studied equation and equations added to it in order to form the automor- 
phic system. In our case we require the compatibility between two equations 
(14. 3 j) which gives restrictions on the function F(t,Ut,p) m the right-hand 
side of the second equation in (|2.4|) . A new feature in our modification of 
the method is that we do this in an explicitly invariant manner by using the 
invariant cross-differentiation 00013110] involving the operators S, Y and 
Y . The resulting five resolving equations have the form 





5(F) = [k(A 


+ A) -5u t ]F 


(5.1) 


F(Y(X) 


-Y(X)) = (r 


+ u tP ){X-X) 


(5.2) 




5(X) = Y{t) 


+ 2u t X - kX 2 


(5.3) 




5{X) = Y{r) 


+ 2u t X - nX 2 


(5.4) 


f(y(X) + 


Y{X)) = -(r 


+ u t p){X + X) 




+ 2k\5(t) 


+ 4u t r + 2F 


+ up 2 + 2u\p 


(5.5) 



The definitions of A, A, r used here are given by equations (|4.4|) . (|3.3|) . 

The resolving equations (|5.1|) . (|5.2j) . ()5.3|) . (|5.4|) . ()5.5|) form a closed re- 
solving system where the second order differential invariant rj = F and the 
third order differential invariants A, A and r are unknown functions of the 
independent variables t,u t ,p. 

An explicit form of the equations (J5.1|) — (J5.5|) is obtained if we use the 
projected operators of invariant differentiation (see (|3.1|) . (|4.3|) . (|4.4|0 

5 = d t + (rep - )d Ut +rd p , Y = d Ut + Xd p , Y = d Ut + Xd p (5.6) 

The commutator relations ()4.2|) were satisfied identically by the operators 
of invariant differentiation. On the contrary, for the projected operators (|5.6j) 
these commutation relations and even the Jacobi identity are satisfied only 
on account of the resolving equations. More than that, the theorem proved 
in jH] claims that commutator algebra (|4.2j) of the operators of invariant 
differentiation 5, Y, Y together with the Jacobi identity gives a commutator 
representation for the resolving system. 
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6 Solutions of the resolving system 
and of the Boyer-Finley equation 

To find particular solutions of the resolving system, we impose various sim- 
plifying constraints. The most obvious ones, like Y = Y or F = 0, lead 
to invariant solutions. These we already know, or can obtain by much sim- 
pler standard methods. The weaker assumption that leads to non-invariant 
solutions is that the operators Y and Y commute 



[Y,Y] = 



T 



-u t p 



(6.1) 



but Y + Y, 



i.e. 



A ^ A and also F ^ 0. With this Ansatz we find the 



particular solution of the resolving system [8] (assuming 2np — uf > 0) 



t = —u t p, A = Ku t + i\J2np 



X = nu t — i\J2np — uj 



(6.2) 



and the expression for F [B] is not needed for obtaining solutions of (jl.lj) . 

To reconstruct solutions of the Boyer-Finley equation starting from the 
particular solution ()6.2|) of the resolving system we use the procedure of in- 
variant integration which amounts to the transformation of equations to the 
form of exact invariant derivative jH] . Then we drop the operator of invariant 
differentiation in such an equation adding the term that is an arbitrary ele- 
ment of the kernel of this operator. This term plays the role of the integration 
constant. This procedure is described in detail in |51I^ITUj. 

Here we present only the final result for solutions of the Boyer-Finley 
equation [H] (see also [Hj in relation to f)6.4[) ) 



u(t, z, z) = In 



(t + b(z))c'(z) 



u(t, z, z) 



In 



c(z) + c(z) 
(t + b{z))c'{z) 



1 + c(z) 



for k = 1 



for k 



(6.3) 



(6.4) 



Here b(z) and c(z) are arbitrary holomorphic functions. One of them is 
fundamental and labels a particular orbit of solutions. The other one is 
induced by a conformal symmetry transformation and can be transformed 
away. These solutions for generic b(z) and c(z) are non-invariant jH]. 
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7 Hodograph transformation of the Boyer- 
Finley equation as a shortcut to its special 
solution 

Hodograph transformation is the interchange of roles of the unknown u and 
one of the independent variables t: u — u(t,z,z) \—*t = t(u,z,z). The 
transformed equation (jl.lj) has the form 

There is an obvious simplifying differential constraint 

tzU - ne u = (7.2) 

valid only for k = 1, since t z t s > 0. Indeed, together with (j7.1|) it implies 
linear dependence of t on z and z and hence t z = e u / 2 + ta ( u ) ) ^_ = e «/2-«a(M) 
with the final result for the solution of (|7.1j) 

t = e u/2+ia{u) z + e u/2 ~ ia{u) z + (7.3) 

where a(u), h(u) are arbitrary functions. This is the solution given in [T7] . 

8 Symmetries of non-invariant solutions 

Point symmetries of the equation (|7.1|) coincide with the symmetries (j2.1j) 
of the original equation (Jl.lj) because the hodograph transformation belongs 
also to a class of point transformations. Since the solution (J7.3|) is the general 
solution of the system of differential equations ([7.1)1 . ([7.2)1 . the symmetries of 
this solution coincide with the point symmetries of that system. A computer 
package 'Liepde' used together with the PDE solver 'Crack' by Thomas Wolf, 
that runs under REDUCE 3.6 or 3.7, gives the symmetry generator of the 
solution ([7.3)1 depending on eight arbitrary functions of u 

X = a(u)d u + [2zg(u) + 2ze u b(u) + te(u) + h(u)]d t (8.1) 
+ [tb(u) + zc(u) + d(u)]d z + {z[2e(u) - a{u) - c(u)] + f(u) + tg(u)e' u }d 2 

Since the solution itself depends only on two arbitrary functions a(u), h(u), it 
is clear that in the set (J8.1)) there exist such symmetries with respect to which 
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the solution (J7.3)) is invariant for any fixed a(u),h(u). Those symmetries 
are not symmetries of the Boyer-Finley equation, therefore (|7.3|) determines 
conditionally invariant (though non-invariant) solution by the definition of 
Levi and Winternitz [TTi] . 

In a similar way we study symmetries of the solutions (j6.3J) and ()6.4|) of 
the Boyer-Finley equation (jl.lj) obtained by Calderbank and Tod ^3] and in 
our paper jHj. They are general solutions of the system of PDEs consisting 
of the original equation and the differential equations 



plus complex conjugate to the last equation, that follow from ()6.2|) by using 
the definitions (|3.3j) . (|2.3j) of r, A, A, p. The symmetries of solutions (|6.3j) 
and (J6.4j) with arbitrary (not fixed!) functions b(z),c(z) coincide with the 
symmetries of this system of four PDEs which turn out to be exactly the 
same as symmetries (j2.1|) of the original equation Hence, our solutions 

with arbitrary fixed generic b(z), c(z) have no symmetries jHJ- Therefore they 
are not conditionally invariant and thus the method of group foliation seems 
to be the only regular method which could give such solutions. 

To prevent misunderstanding, we warn against mixing up the invariance 
of a particular solution with a fixed choice of arbitrary functions with the 
invariance of its orbit when arbitrary functions are allowed to transform under 
a symmetry transformation since the orbit is always an invariant manifold. 

9 Transformations of differential invariants 
and invariant differential operators 

Our main problem is to understand how we could arrive at the solution (|7.3|) 
by the method of group foliation, i. e. to determine the orbit of this solution. 

The first step is to make the hodograph transformation of the whole group 
foliation structure. In this section we make the hodograph transformation of 
differential invariants and of invariant differential operators and in the next 
section we present the hodograph transformation of the resolving equations. 

Since the hodograph transformation is a point transformation, it con- 
serves point symmetries, differential invariants and operators of invariant 
differentiation. We change the definitions of differential invariants p, 77 and 



U Z zt — 0, 



u Z zz = u z u Z z + e u u zt KU t + iJ 2ne 
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invariant differential operators after the hodograph transformation and there- 
fore from now on we label the old invariants and operators with the tilde. 

The set of hodograph-transformed differential invariants of up to the sec- 
ond order inclusive consists of the new unknown t = t(u, z, z) and its deriva- 
tives t u = 1/ut, t uu = —Uft/uf and two more invariants 

P 6 \pztz~tuu ~\~ t'upzz t u (tzt<u z -\- tzt-uz)\; P ^ ^22 P/^u 
T) ^ (t%tuu t u t uz ^j {tzt uu lid uz): V ^ ^ztUzt V/^u 

The explicit invariant form of the transformed Boyer-Finley equation (|7.1|) 
becomes 

tuu = Kp + t u (9.1) 

Let us denote by D u the total derivative with respect to u taken for the 
constant values of z, z. Since the invariant differential operator 6 — D t — 
(l/t u )D u is a total derivative with respect to invariant t, its definition is not 
changed: 5 = 5. The definitions of A, A become (A = — (l/tfyA, A = 

-(VOA) 

We choose t,t u ,p as the independent invariant variables in accordance 
with our previous choice. The action of new operators of invariant differen- 
tiation on the independent variables becomes 

5(t) = 1, 5{t u ) = k£ + 1, A(t) = A(t) = 0, A(t u ) = A(t u ) = V 

5(p) = r, A(p) = a, A(p) = o (9.2) 

where for 5{t u ) we have used the equation (|9.1jl . 

New invariant differential operators Y = (l/i])A and Y = (l/r])A have 
the properties 

Y(t) = 0, Y(t)=0, Y(t u )=Y(t u ) = l, Y(p) = X, Y(p) = \ (9.3) 



10 Hodograph transformation of resolving 
equations 

We consider the automorphic equation 

V = A(t u )=F(t,t u ,p) (10.1) 
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together with the transformed Boyer-Finley equation (|9.1|) and derive the 
resolving equations as compatibility conditions of these two equations. The 
resulting five equations coincide with the hodograph transform of the old re- 
solving equations with the only change caused by the difference in definitions 
of differential invariants and invariant differential operators 



5{F) = -}-[k(\ + \) + 1)F 
F{Y{\)-Y{\)) = {\-\) 
5(X) = Y(t) + k>j- (j- — A 



5(\) = Y(r) + K-(^-X 
F(F(A)+F(A)) = (A + A) 
+ 2 K t 2 u 5(r) - (12p + 4/rf u ) 



P 2 2 „ ^ 
t u T-3K^- + -F-2p 



(10.2) 



F p 2 
-t u r + 4 — + 3k — + 2p 

tii tu j 



T + 



tl 



K 



P__P_ 

tl u 



Thus, we have 4 unknowns F, A, A and r and 3 independent variables t, u t , p. 
Here 8,Y,Y denote the projected operators (see (|9.2|) . ()9.3|) ) 

8 = d t +^K^ + l^d tu +rd p , Y = d tu + Xd p , Y = d tu + Xd p (10.3) 



11 Towards the invariant description of the 
hodograph solution 

Our final goal is to obtain a complete set of independent additional relations 
between differential invariants which are satisfied on the hodograph solution 
of ^7j. Combining them with the resolving equations we would arrive at 
the particular solution of the resolving system which selects the orbit of this 
solution and makes it possible to reconstruct the solution ()7.3|) . Using these 
invariant relations with the opposite sign of k, we hope to obtain new non- 
invariant solutions of the Euclidean Boyer-Finley equation (Jl.lj) with k — — 1 
which are more physically interesting. 
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Up to now we have succeeded to discover only three additional relations 
between differential invariants satisfied identically on the solution (|7.3|) 




u 



4 



1 



[(2 P + t u )(\ + \)-iyj4$-tl 



(A - A)] 



(11.1) 



A($) 



4 



i 



v /4$-t2| 2 [$(A + A) + (p 2 + t uP )(A 



A) 



t 2 u + i{2p + t u )y/4$-tl } A) + f $F 



(11.2) 



together with the complex conjugate to (jll.2j) where $ = F — p 2 — t u p. We 
use the first one of these equations to exclude r from the remaining equations. 
It turns out that these additional equations are not sufficient for solving the 
resolving equations (jl(J.2j) . We are now in the process of searching for missing 
relations. 



In the first part of this paper we have given a short exposition of our de- 
velopment of the method of group foliation as a method for obtaining non- 
invariant solutions of PDEs which admit an infinite dimensional symmetry 
group. The method was illustrated by the physically interesting Boyer-Finley 
(or heavenly) equation. It was shown how the constraint of commutativity 
of two operators of invariant differentiation led to non-invariant solutions of 
the Boyer-Finley equation. 

In the second part of the paper we discussed the hodograph solution of 
Manas and Alonso. We have shown how it naturally arises after making the 
hodograph transformation of the Boyer-Finley equation without any refer- 
ence to the hydrodynamic-type systems. It turns out that this solution has 
extra symmetries as compared to the equation itself and that, as a conse- 
quence, it is conditionally invariant, unlike non- invariant solutions obtained 
previously. 

Our final goal in this part of the paper was to obtain a complete set 
of independent additional relations between differential invariants which are 
valid for the hodograph solution of and using them to obtain a particular 
solution of the resolving system which selects the orbit of the hodograph 
solution. In this way we hoped to obtain more general solutions of resolving 



12 Conclusions and outlook 
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equations if we would skip one of the additional relations between differential 
invariants and, as a consequence, more general solutions of the Boyer-Finley 
equation in a hodograph form. More than that, those extra relations should 
not be sensitive to the sign of k as it was in the case of solutions ()6.3|) and 
()6.4p and hence we could reconstruct solutions of the equation (|7.1|) with 
the opposite sign k — — 1 which was impossible using the method of the 
paper [T2j. Until now we have not yet succeeded to fulfil this program and 
present here only partial results. 

Acknowledgements. I thank Evgenii Ferapontov and Yavuz Nutku for 
useful discussions. 

References 

[1] L.V. Ovsiannikov. Gruppovoi analiz differentsial'nykh uravnenii. M.: 
Nauka, 1978. English translation: Group Analysis of Differential Equa- 
tions. New York: Academic, 1982. 

[2] P. Olver. Applications of Lie Groups to Differential Equations. New 
York: Springer- Verlag, 1986. 

[3] P. Winternitz. Group theory and exact solutions of nonlinear partial 
differential equations. In: Integrable Systems, Quantum Groups and 
Quantum Field Theories. Dordrecht: Kluwer, 1993. P. 429-95. 

[4] S. Lie. Uber Differentialinvarianten. Math. Ann. 1884. V. 24. P. 52-89. 

[5] E. Vessiot. Sur l'integration des sistem differentiels qui admittent des 
groupes continus de transformations. Acta Math. 1904. V. 28. P. 307- 
349. 

[6] Y. Nutku, M.B. Sheftel. Differential invariants and group foliation for 
the complex Monge-Ampere equation. J. Phys. A: Math. Gen. 2001. V. 
34. P. 137-156. 

[7] Y. Nutku, M.B. Sheftel. Group foliation approach to the complex 
Monge-Ampere equation. Teoreticheskaya i Matematicheskaya Fizika. 
2001. V. 127, No. 3. P. 465-474 (in Russian). English translation: Theor. 
Math. Phys. 2001. V. 127, No. 3. P. 808-816. 



14 



L. Martina, M.B. Sheftel, P. Winternitz. Group foliation and non- 
invariant solutions of the heavenly equation. J. Phys. A: Math. Gen. 
2001. V. 34. P. 9243-9263. 

M.B. Sheftel. Method of group foliation and non-invariant solutions of 
invariant equations. Proceedings of Institute of Mathematics of NAS of 
Ukraine. 2002. V. 43, Part 1. P. 215-224. 

M.B. Sheftel. Method of group foliation and non-invariant solutions of 
partial differential equations. Example: the heavenly equation. Eur. 
Phys. J. B. 2002. V. 29. P. 203-206. 

CP. Boyer, J.D. Finley III. Killing vectors in self-dual, Euclidean Ein- 
stein spaces. J. Math. Phys. 1982. V. 23. P. 1126-1130. 

I.A.B. Strachan. The dispersive self-dual Einstein equations and the 
Toda lattice. J. Phys. A: Math. Gen. 1996. V. 29. P. 6117-6124. 

T. Eguchi, P.B. Gilkey, A.J. Hanson. Gravitation, gauge theory and 
differential geometry. Phys. Rep. 1980. V. 66. P. 213-393. 

D. M.J. Calderbank, P. Tod. Einstein metrics, hyper complex structures 
and the Toda field equation. Differ. Geom. Appl. 2001. V. 14. P. 199-208. 

Y. Nutku, M.B. Sheftel. A family of heavenly metrics. 2002. Preprint of 
arXiv.org e-Print archive arXiv:gr-qc/0105088, 

R.S. Ward. Einstein- Weil spaces and SU (oo) Toda fields. Class. Quan- 
tum Grav. 1990. V. 7. P. L95-L98. 

M. Manas, L.M. Alonso. A hodograph transformation which applies 
to the heavenly equation. 2002. Preprint of arXiv.org e-Print archive 
|arXiv:nlin.S I /0209050) 

E. V. Ferapontov, D.A. Korotkin, V.A. Shramchenko. Boyer-Finley 
equation and systems of hydrodynamic type. Clas. Quantum Grav. 2002. 
V. 19. P. L1-L6. 

D. Levi, P. Winternitz. Non-classical symmetry reduction: example of 
the Boussinesq equation. J. Phys. A: Math. Gen. 1989. V. 22. P. 2915- 
2924. 



15 



[20] A. Tresse. Sur les invariants differentiels des groupes continus de trans- 
formations. Acta Math. 1894. V. 18. P. 1-88. 



16 



